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ABSTRACT 



Aims. In this paper, we investigate the behaviour of the number of Lyman continuum ionizing photons as compared to the actual number of 
hydrogen recombinations in HII galaxies. 

Methods. We evaluate the number of ionizing photons from the population synthesis of spectra observed in the visible, extrapolating the spectra 
to the extreme ultraviolet (EUV), beyond the Lyman limit. We check for possible systematic deviations of the predicted ionizing spectra in the 
EUV by comparing the ratio of the predicted number of ionizing photons to the number of recombinations, as measured in H/?, A log Q(H°), 
with the metallicity. We find that, as far as the number of ionizing photons is concerned, no systematic tendency can be detected. 
The H/J equivalent width can be understood as a nebular age indicator, decreasing with age, although the observed H/3 equivalent width can 
also be affected by the contribution to the continuum by the accumulation of previous, non-ionizing stellar populations. 

Results. We attribute the increase of A log Q(H°) with the age of the burst to the fact that more and more ionizing photons escape the nebulae 
when the nebulae get older, because of their increasing, expansion-induced subfragmentation. 



Key words. Galaxies: starburst - Galaxies: ISM - ISM: HII regions - Line: formation 



1. Introduction 

In the giant HII regions present in HII galaxies, photoionization 
converts ionizing photons of the Lyman continuum in the ex- 
treme ultra-violet (EUV) (e.g., with A < 912A or E > 13.6 eV) 
into emission lines, ranging from the ultra-violet to the far in- 
frared; the most prominent ones are recombination lines and 
collisional lines of atoms, generally in their first and second 
ionization stages. For photons in the EUV, the radiative world 
is very different from the one for photons with higher wave- 
lengths. In the latter, small amounts of neutral hydrogen al- 
ready cause high opacity, and photons are rapidly absorbed. 

The EUV part of the ionizing stars in HII regions is poorly 
known; it is inaccessible observationally and can only be 
known from stellar atmosphere models. Ionizing stars in HII 
regions, mainly O and B stars and sometimes Wolf-Rayet stars, 
possess strong winds, and strong non-ETL effects are pr esent in 
their atmospheres (Hillie r & Miller £l998); Pauldrach bOOll) : 
Schaerer & Vacca dl998h : Schaerer & de Koter J 19971) : Smith et 
al. (120021) '). In the EUV, various bound-free ionization disconti- 
nuities are present, strongly affecting the spectra; however, dif- 
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ferences b etween models by different authors can be significant 
(Morisset(2003)). 

The main constraint on these model spectra comes from 
the reconversion of their energy in to em ission lines, which 
are more easily observed. Morisset (2004), for instance, used 
high ionization stages lines observable in the far infrared 
to try to constrain the spectral energy distribution of ioniz- 
ing stars of HII regions in the EUV. This kind of work re- 
quires the use of photoionization codes , like CLOUDY (e.g. 
Ferland ill 9861)1 NEBU (e.g. Pequignot jl980V ). PHOTO (e.g. 
StasinskaUSi). 

However, more basic information can be extracted from 
just the analysis of hydrogen recombination lines: in the pho- 
toionization equilibrium, the number of ionizations of neutral 
hydrogen is balanced by the number of recombinations. The 
number of ionizations can easily be calculated from the (theo- 
retical) ionizing spectrum, and the number of recombinations 
can be counted in hydrogen recombination lines that hardly 
depend on the plasma parameters (e.g. electronic temperature). 
The intensity of hydrogen recombination lines should thus 
simply be proportional to the integrated luminosity of the 
ionizing spectra in the EUV. This has been known since the 
pioneering works by Zanstra Jl93lh on planetary nebulae, but 
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the measurements of EUV luminosities by hydrogen Balmer 
recombination lines strongly rely on the hypothesis of the 
optical thickness of the nebulae to Lyman continuum photons. 
In planetary nebulae, Tylenda & Stasinska (1989) have shown 
that the optical thickness strongly affects the measurements of 
their EUV luminosities. 

As far as HII regions are concerned, direct measurements 
of the Lyman continuum for local starburst galaxies show that 
either they are optically thick or only a low fraction (< 1 0%) o f 
ionizing photon s are a llowed to escape ( Leither er et al. ([1995); 
Heckm an et al. d200ll) : Deharveng et al. J2001I) : Bergvall et al. 
(2006)). On the other hand, also from direct Lyman continuum 
measurements, Steidel et al. d200ll) found no evidence for the 
absorption of ionizing photons in high redshift (z^3.4) Lyman 
Break Galaxies; they estimate a lower limit for their Lyman 
continuum escaping photons of 50%. 

In the local universe, indirect measurements of lower lim- 
its of the escaping fraction of ionizing photons are possible 
from the emission of the diffuse ionized gas (DIG): Zurita et al. 
(2000) estimate that the luminosity of the DIG in local spiral 
galaxies can reach 25-70% of their total Ha luminosity. One 
has to bear in mind that the actual escaping fraction can be 
higher, as some photons may leak sideways, perpendicularly 
to the gal actic p lane. In the Milky Way, Bland-Hawthorn & 
Malloney (119991) estimate the upper bound to the leaking frac- 
tion of ionizing photons to be only of the order of 10%, from 
the ionization of clumps in the magellanic stream. 

The understanding of the emission lines sequence in HII 
galaxi es can bring some insight into this problem; Stasinska 
et al. l l200ll) and Stasinska & Izotov (2003) concluded from 
the statistical photoionization modelling of large samples (ei- 
ther published catalogues from spectrophotometric surveys or 
SDSS early data release ') that some systematic effect was 
present in their optical thickness, e.g., that a fraction of HII 
galaxies was optically thin, and suggested that HII g alaxies 
were becoming thinner with time. Castellanos et al. (2002) 
reached similar conclusions from the detailed photoionization 
analysis of 3 giant HII regions in nearby galaxies, e.g., that 
they were losing a significant fraction of their ionizing photons 
(between 10-70 %). 

Photoionization modelling, however, generally does not 
take into account the presence of dust, which can absorb 
Lyman continuum photons before they even have a chance 
to ionize the gas and then re-emits them in the infrared. A 
significant fraction of the photons that are interpreted as 
escaping photons may actually be used to heat up th e dust . 
Assuming that no ionizing photon escapes, Inoue ll200ll) 
infers from the comparison of radio (gas free-free emission) to 
infrared (dust emission) luminosities of individual HII regions 
that the fraction f, of ionizing photons used for photoionization 
varies between 0.4 and 0.7. Since he assumes that all ionizing 
photons are used either to ionize the gas or heat up the dust, 
the fraction 1 - f; used to heat up the dust should only be 
considered as an upper limit, as some ionizing photons may 
simply escape. Nonetheless, this upper limit, of the order 

1 Sloan Digital Sky Survey (www.sdss.org) 



of 0.6, is significant. Inoue et al. (2001) and Hirashita et al. 
(2003) confirmed this result for spiral galaxies, finding upper 
bounds for 1 - /, ranging from 0.4-0.7 for 1 - ff. 

It has to be emphasized that the problem of the optical 
thickness of the ionized gas in HII galaxies is not restricted 
to these galaxies. Reionization by the first hypermassive pop- 
ulation III stars is a connected question, as the gas ionization 
responds to the energy input and to the incoming radiation 
in nearly the same way as in local HII galaxies. Most likely, 
the gas constituting the intergalactic medium at this stage of 
the evolution of the Universe is heterogeneous and clumpy. 
Its interactions with primordial stars can only be understood 
in terms of complex 3D hydro dynamical and photoionization 
modelling (Abel et al (2000)), whereas observational con- 
straints are nearly non-existant. HII galaxies might be the 
closest that we have to a laboratory of the physical conditions 
of reionization. 

Effects originating from the optical thickness, or geomet- 
rical effects, are intrinsically hard to distinguish from effects 
arising from uncertainties concerning the ionizing spectra or 
from the presence of dust. One way to examine the question is 
to look at the metallicity of the ionizing clusters. The integrated 
spectra of an ionizing population should vary with metallicity, 
partly because of the variation of the metal opacities of the 
stars that produce the ionizing spectra, but mostly because 
of the variations of evo lutionary tracks due to the metallicity 
(Schal ler et al dl992l) : Charbonnel et al. dl993l) : Schaerer 
et al. dl993al) . dl993bl): Mow l avi et a l. dl998l) ; Girar di et a l. 
dl996l) : Fagotto et al. Jl994ah . |l994b); Bressan et al. dl993l) ). 
A great part of the opacity in the atmospheres of hot stars 
in the EUV is due to bound-free transitions of hydrogen and 
heliu m - and thus is not affected by metallicity (see Leitherer 
et al. dl999» ). The amount of dust i s also expected to increase 
with metallicity (Inoue et al d200lh . Variations of the optical 
thickness with metallicity (which can easily be measured from 
oxygen emission lines of the ionized gas) should therefore be 
indicative of systematic effects of the ionizing spectra, or of 
the dust, and represent an important check on how to segregate 
effects of optical thickness from other ones. HII galaxies are 
actually quite fit for this purpose, because their HII regions 
are giant, with hundreds or thousands of ionizing stars; in this 
case, it is much easier to predict the ionizing radiation with a 
population synthesis code than in normal HII regions, which 
only possess a few ionizing stars, because the counting noise 
strongly affec ts the predicted ionizing spectrum (Cervino & 
Valls-Gabaud d2003ll ). 

In this paper, we investigate effects of optical thickness in 
a sample of HII galaxies by comparing the number of ionizing 
photons produced by stars and the number of recombinations. 
In Sect. 2, we will describe the sample, and our method to eval- 
uate the number of ionizing photons from the stellar spectra. In 
Sect. 3, we will describe how we can understand statistically 
the evolution of the (mean) optical thickness with time, e.g., 
with the nebular age, and how we can check for systematic 
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effects of the model ionizing spectra with metallicity. Finally, 
Sect. 4 will present the results and a summary. 

2. Sample and methodology 

Our study is based on the spectro photometric catalogue of HII 
galaxies by Kehrig et al. (2004). This catalogue contains ho- 
mogeneous longslit observations of 188 spectra of individual 
regions in a sample of 1 1 1 HII galaxies. The choice of the re- 
gions was made by luminosity, the slit being positioned along 
the brightest regions. These regions are typically made up of 
bright knots, generally encompassing several stellar superclus- 
ters, with sizes of a few hundreds of pc. 

The Kehrig et al. catalogue includes most of the southern 
hemisphere HII galaxies listed in Terlevich et al. J 1 99 lit - but 
with a higher signal/noise. It also includes nearby HII galaxies 
that have been discovered since the publication of Terlevich et 
al.'s catalogue. The stellar population content of the galaxies 
from this catalogue, as derived from their spectra, has already 
been discussed in Wester a et al. ll2004l) . 

Westera et al. il2004l) based their study on the analysis of 
modified Lick indices (Worthey ( 1994); Worthey et al. ( Il994l) '). 
What we want to do here is to use a similar population synthe- 
sis method, based on a full spectra fitting in the 4000-7000 A 
range, with the purpose of predicting the ionizing spectrum in 
theEUV, for A < 912A. 

A best-fit procedure was employed using three different li- 
braries of simple stellar populations (SSPs). The first SSP li- 
brary (hereafter the "BC99" library) was produced using the 
Bruzual and Chariot 2000 Galaxy Isochrone Spectral Sy nthesi s 
Evolution Library ( GISSE L ) code (Chariot & Bruzual Jl99ll) : 
Bruzual & Chariot dl993l) . i2000lf>. imp lementing the Padova 
2000 isochrones (Girardi et al. (2000)) combined w ith th e 
BaSeL 3 1 "Pad ova 2000" stellar library (Westera et al. J200ll) : 
Westera (2001)). The second SSP library, "Starburst", consists 
of sp ectra f rom the STARBURST99 data package (Leitherer 
et al. dl999l) '). using the option of including nebular continuum 
emission (Fig. 1 on the STARBURST99 we b site) It imple- 
ments the BaSeL 2.2 library (Lejeune et al. ill997ll . £998)), 
and for stars with strong mass loss, it also take s into account ex- 
tended model atmospheres by Schmutz et al.Jl992), combined 
with the Geneva isochrone s (Mey n et et al. ( 119941) : Schaller et 
al. Jl992h : Schaerer et al. Jl993al) . Jl993bl) : Charbonnel et al. 
( 1993J)). For old populations, the "BC99" spectrum was used, 
since the Starburst99 data package only contains spectra up to 
900 Myr. In addition t o thes e two libraries - which were already 
used in Westera et al. ( 2004) - we also used a library with higher 
spectral resolution, "BC03", produced by employi ng the 2003 
version of the GISSEL code (Bruzual & Cha riot J2003I) ') and 
the Padova 1995 isochrones (Fagotto et al. ( 1994b); Girardi 
et al. i 19961) ) combined with the STELIB (Le Borgne et al. 
J2003I) ') stellar library. As the nebular continuum emission is 
important for the spectral shape of young populations, it was 
also added to the spectra in the "BC99" and " BC03" libraries, 
in the way described by Leitherer et al. (1999). 
Table^provides a summary of the evolutionary tracks and stel- 
lar libraries used for the three SSP libraries. 



Table 1. Evolutionary tracks and stellar libraries used in the 
SSP libraries. 



SSP library 


Tracks 


Stellar library 


"BC99" 


Padova 2000 


BaSeL 3.1 "Padova 2000" 


"Starburst" 


Geneva 


BaSeL 2.2 


"BC03" 


Padova 1995 


STELIB 



Table 2. Wavelength ranges that were not used for the spectral 
fit. 



Range (A) 


"Contamination" source 


3960-3975 


He+[NeIII]3967 


4095.25-4106.5 


H<5 


4331-4370 


Hy+[OIII]4363 


4465-4479 


HeI4471 


4848-4871 


H/J 


4943-5020 


[OIIIJ4959+5007 


5446-5574 


telluric lines 


5867-6320 


HeI5876, [OI]6300, [SIII]6312 


6375-6470 


telluric lines 


6520-6600 


Ha+[NII] 


6700-6740 


[SII] 6717, 6731 



Since for the present purpose, it is more important to re- 
produce the shapes of the galaxy spectra than to find the actual 
parameters of the sub-populations, we fit the full spectra from 
3910 A to 6880 A instead of single spectral indices, except 
for those various parts of the spectra showing "contaminations" 
from different emission and/or telluric line sources. Table|2]lists 
all the regions that were not used for the fits. 

In the present paper, we modelled the actual population as 
being composed of an old, an intermediate, and a young stellar 
population (in Westera et al. (2004), we only decomposed the 
actual population as a young +intermediate one and an old one). 
The characteristics and free parameters of the three populations 
are summarised in Table [3] 

The spectra were first blueshifted to their rest wavelength. 
Internal gas extinction was then corrected for, using the ob- 
served H a /Hp emission lines ratio. The resolution of our spec- 
tra was sufficient (^ 5A) to easily deconvolve Hq, f rom th e 
[MI] doublet at ,L16548, 6583A (see Kehrig et al. ( 12004 1. 
As explained in next section, we did no t adopt the emission 
line intensities from Kehrig et al. II2004I) . but we remeasured 
them after fitting the stellar continuum by population synthe- 
sis. In HII galaxies spectra, H a and H^ are present both in 
absorption (from the stars) and in emission (from the ionized 
gas); our procedure allowed us to disentangle both compo- 
nents. The repartition of gas and stars inside a galaxy lets them 
be affected in different manners by extinction (see Calzetti et 
al. (11994)). We dereddened the stellar continuum spectra em- 
ploying the formulas E(B - V) = 0.665 x Chp (Sampson 
(2003)) and E(B - V) co „, = 0.44 x E(B - V) gas (thus, 
E(B - V) cont = 0.2926 x C H p) and the extinction law of Fluks 
Jl994l) . The factor 0.44 was introduced to correct for system- 
atic differences between the extinction as deriv ed from the stars 
and from the ionized gas (Calzetti et al. (119941) 1. 
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Table 3. Possible values of the population parameters. 



Parameter 


Possible values 


(M y + Mj):M 


0:1, 1:100, 1:30, 1:10, 1:3, 1:1, 3:1, 10:1, 1:0 




0:1, 1:30, 1:10, 1:3, 1:1, 3:1, 10:1, 30:1, 1:0 


age y 


1, 2, 3, 4, 5, 6, 7, 8, 9 Myr 


age/ 


10, 20, 50, 100, 200, 500 Myr 


age„ 


fixed at 5 Gyr 


[Fe/H] v = [Fe/H], 


taken from Cuisinier et al. (2006) or -1 


[Fe/H] 


fixed at -1.5 



The best- fitting population was found by a% 2 algorithm. To 
be able to calculate the x 1 estimator, the empirical and theoret- 
ical spectra have to be on the same wavelength grid (at least 
in the range used for the fit). This was done by rebinning the 
empirical spectra to the resolution of the theoretical spectra us- 
ing a Gaussian kernel function with a full width at half max- 
imum (FWHM), corresponding to the resolution of the theo- 
retical spectra (20 A for "BC99" and "Starburst" and 1 A for 
"BC03"). 

In the end, a selection of acceptable solutions was made. All 
spectra that had too bad a signal-to-noise or a strange shape, in- 
dicating calibration problems in the data reduction, were elim- 
inated from the sample. This selection was done by eye. It left 
us with solutions for 105 spectra of individual regions in 72 
galaxies. 

The best-fitting population parameters for the fits with the 
"BC03" library can be found in columns 4 to 8 of Table IA.fi 
Figure ^ shows an example of an observed spectrum and the 
resulting composite spectra in the wavelength range where the 
fits were made, and Fig. [2] represents the extrapolations to the 
UV and EUV. 

The introduction of a third partial population again serves the 
purpose of better reproducing the spectral shapes. In reality, 
the quality of our spectra is not high enough to determine 
the parameters (metallicities, ages, and masses) of three sub- 
populations, so the parameter values given in Table lA.il should 
be taken with a grain of salt. However, the parameters tha t can 
directly be compared with those derived in Westera et al. ( 2004) 
agree reasonably well: We obtain the same value for t he par am- 
eter (M y + Mi):M as for (M y+i ):M in Westera et al. ( 2004) for 
38% or 49% of the spectra (depending on the spectral library, 
see next paragraph), and in another 36% or 28%, these two pa- 
rameters differ only by one "step" in the grid. In 49% or 46% 
the luminosity-weighted (geometrical) mean of age y and age, 
differs by less than two "steps" from the age y+ i calculated in 
Westera et al. (2004), which corresponds to a factor of less than 
5 in age y+i . 

In Westera et al. J2004I) . we discussed the diffe rence s 
arising from adopting the Starburst 99 (Leitherer et al. (1999)) 
or the BC93 (Bruzual & Chariot i 1 9931) 1 codes for the young 
population (the old one being modelled with the BC93 code 
in both cases). In this paper we enhance these results by 
modelling both the old and young populations with the BC03 
(Bruzual & Chariot (2003)) version of their code. Each 
of these codes, or code versions, has their advantages and 



inconveniences: Starburst 99 is particularly fit for the study 
of young populations in the visible/ultraviolet, but does not 
extend to ages older than lGyr; BC93 is made to model all 
kinds of populations, but is optimized for old ones. Finally , 
the BC03 version of the code (Bruzual & Chariot (2003)) 
includes empirical spect ral libraries of higher resol ution, 
BC-2000 (Pickles Jl998h : Bruzu al & Chariot ll2003l) 1 and 
STELIB (le Borgne et al. lT2003h . and, therefore, is better 
suited for calculating Lick indices (because of its higher 
spectral resolution) than Starburst 99 or BC93. However, it 
has a poor coverage of parameter space for the ionizing hot 
stars. Furthermore, the extension of the BC03 version to wide 
and comprehensive wavelength ranges is made by combining 
STELIB spectra in the visible, Pickles library ( i 19981) 1 in the 
UV and IR, and BaSeL theoretical spectra beyond 2.5 fim of 
for wavelengths shorter than 1200 A. Comparing the results of 
different codes or versions, we hope to be able to better assess 
possible systematic effects. 



3. Evaluation of the photon loss rate 

We evaluate the number of ionizing photons with the method 
exposed in the previous section from an extrapolation of the 
ionizing spectrum to the EUV, based on population synthesis. 
The number of recombinations can also be calculated from a 
recombination line of hydrogen, like H/3. If the number of pre- 
dicted ionizing photons is correct, in the case of an optically 
thick ionized nebula, both numbers are equal (in the photoion- 
ization equilibrium). In an optically thin nebula, however, not 
all the ionizing photons will actually be able to ionize the gas, 
and some will escape. The ratio of the number of (predicted) 
ionizing photons to the number of recombinations can thus be 
considered an indicator of optical thickness. We define: 



AlogQ(H°) = logQ(H°) ion - logQ(H ) I1 



(1) 



where Q(H°)i on is the number of ionizing photons, as deter- 
mined from the extrapolation of the ionizing spectrum to the 
EUV, based on population synthesis, and Q(H°) rec is the num- 
ber of recombinations as calculated from the H/3 recombination 
line, after correction from absorption, as estimated from Ch/s- 
To properly measure the equivalent width of the H/3 emission 
line, we had to correct its strength for the reduction it suffers 
from the underlying absorption line. For this purpose, we first 
made a best fit of the continuum in the 3910-6880 A wave- 
length range, in a similar way to Westera at al. (2004), but 
putting extra weight on the H6 index (the wings of the H6 ab- 
sorption line) and using the theoretical spectra from the BC- 
2000 stellar library. Although this library only exists for solar 
metallicity, its resolution is much higher than for the other li- 
braries, and the best-fitting spectra reproduce the shape of the 
absorption lines better than the fits obtained with the other li- 
braries. Note that these fits were only made for calculating the 
H/3 equivalent widths and have nothing to do with the popu- 
lation synthesis described in Sect. 2. After subtracting the (re- 
binned) best-fitting spectra from the galaxy spectra, we could 
measure the strength of the emission line and thus calculate the 
equivalent widths EW(Hp). This is illustrated in Fig. [3] Results 
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Fig. 1. Example of a "best fit". The thin (red) lines represent the empirical spectrum (UM 137(W), taken on August, 18, 1998); 
the thicker black lines show the best-fitting spectra obtained using the different SSP libraries. 



are given in column 2 of Table lA. 11 

Using these definitions, AlogQ(H°) will be equal to for op- 
tically thick nebulae and increase for nebulae getting optically 
thinner and thinner. The A log Q(H°) values of individual spec- 
tra can be found in column 3 of Table lA. fl 



As we are more interested here in giant HII regions than in 
the HII galaxies themselves, we will focus our discussi on on 
the spectra of the individual regions from Kehrig et al. 's ll2004 
catalogue, and not on the spectra of the individual galaxies. As 
a matter of fact, we checked if our results were the same con- 
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4850 4900 4850 4900 

A [A] 

Fig. 3. Illustration of how the H/3 emission line equivalent 
width was measured through the example of UM137(CentE). 
In the left panel, the solid line represents the empirical spec- 
trum, whereas the short-dashed line shows the best fit using 
the BC-2000 stellar library. The shaded region between these 
two lines (shaded with lines running from the upper left to the 
lower right) shows the area used to calculate the emission line 
strength, and the shaded region (lines running from the upper 
right to the lower left) between the solid line and zero level (the 
long-dashed line) shows the location of the continuum band. 
The right panel shows the emission line after subtracting the 
(rebinned) best fit. The shaded area corresponds to the shaded 
area between the spectrum and the best fit from the left panel. 
Here, too, the zero level is shown as a long-dashed line to show 
the quality of the subtraction. 

sidering only the integrated spectra of galaxies, and we found 
excellent agreement. 

We could have some systematical bias due to the fact 
that we observed the galaxies, or giant HII regions, through 
apertures defined by the slit, and thus not encompassing the 
HII regions (galaxies) entirely. As we centred the slit on the 
brightest regions, we estimate that the quantity of nebular Hy8 
photons we lose should not exceed a factor of ^ 3 (0.5 in 
logarithm). 



3.1. Evaluation of the ionizing spectra 

It is intrinsically hard to disentangle effects originating from the 
ionizing spectra from effects arising from the optical thickness 
of the nebulae: For a given HII region, an erroneous modelling 
of the number of ionizing photons can give results equivalent 
to the ones that would be attributed to an optically thin nebula. 
The accuracy of the modelling of the stellar spectra in the EUV 



can be assessed by the comparison of line intensities of ion s 
of different ionization potentials, as in Morisset et al. (2004). 
Much more basic information can, however, be extracted from 
systematical tendencies of A log Q(H°) with metallicity. If the 
predicted ionizing spectra are wrong for some reason, their sys- 
tematical excess or deficiency in ionizing photons will vary 
with metallicity because of the intrinsic change of the inte- 
grated ionizing spectra (see Leitherer et al. ( 1999)) and will 
be seen in A log Q(H°), if plotted as a function of metallicity. 
On the other hand, geometrical effects (optical thickness) are 
not expected to vary with metallicity and will thus just produce 
a systematical offset in AlogQ(H°). 

Figure 0] shows a comparison of AlogQ(H°) with the 
oxygen abundance as derived from emission lines of the 
gas (Cuisinier et al. (2006)). Most giant HII regions have 
A log Q(H°) > 0, indicating that they let ionizing photons es- 
cape, since they are optically thin. A huge scatter exists, how- 
ever, with some giant HII regions reaching negative values of 
AlogQ(H°). Nonetheless, as we will show in Sect. 3.3, we be- 
lieve it to be due to uncertainties in the population synthesis fit- 
ting procedure. If we look at the bulk of the giant HII regions, 
however, they lie at a constant value of A log Q(H°), around 0.8, 
with no clear tendency for oxygen abundance. 

Another point that could affect our evaluation of the number 
of ionizing photons is dust absorption. As stated in Sect. 2, we 
took great care in correcting the observed spectra for absorp- 
tion before performing any fitting procedure. We did this by 
measuring the reddening from the observed Ha/H/? line ratio. 
Absorption for the stellar population, however, is not necessar- 
ily the same as for the gas phase: the geometrical repartition is 
not identical in respect to the dust. We applied a correction fac- 
tor on the reddening for th e stella r population to compensate for 
this effect (Calzetti et al. but it is important to check 

that no remaining differential effect exists. Furthermore, dust 
can absorb high energy ionizing photons, affecting our evalua- 
tion of A log Q(H°). 

The dust quantity should increase with metallicity; any ef- 
fect due to an erroneous evaluation of the absorption should 
affect the variation of AlogQ(H°) with metallicity, but, as can 
be seen in Fig. |4] no clear tendency can be detected. 

This might seem contradictory to the results by Inoue et 
al. J200lh : From the analysis of the infrared emission of HII 
regions, due to the dust, and their radio emission, due to the 
gas, they find that the fraction of Lyman continuum photons 
absorbed by the dust increases with metallicity (assuming that 
no photons are leaking). The metallicity range of their HII re- 
gions is different, however, and extends to higher metallicities 
than our HII galaxies, though some overlap exists. A closer in- 
spection of their data shows that at low metallicities, equivalent 
to the ones we encounter in our sample ([0/H]< -0.5), the frac- 
tion of Lyman continuum photons used for dust heating I — ft 
is rather constant at ^ 0.2. It only increases for higher metallic- 
ities. One must bear in mind that, as stated in the introduction, 
since they do not consider the possibility of photons leaking 
out, 1 - fi is an upper limit to the fraction of Lyman continuum 
photons used to heat up the dust. 

Thus, from our analysis of the variation of A log Q(H°) with 
metallicity, (i) we cannot put into evidence any deficiency from 
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the theoretical ionizing spectra, and (ii) we cannot rule out ef- 
fects of absorption by the dust, if they are constant with metal- 
licity. Inoue et al. showed that such an effect exists, 

but at the metallicities of our HII galaxies sample, it should 
be lower than 20%. 

As a consistency check, in Fig.|5]we compare AlogQ(H°) 
with the measured Cup values; effects due to the dust, provok- 
ing either an erroneous evaluation of the absorption for the stel- 
lar population in the visible or an absorption of high energy 
ionizing photons, should vary with the dust quantity, e.g., the 
reddening. As can be seen in Fig. [5] no tendency can be de- 
tected. 

3.2. Temporal evolution of the optical thickness to the 
ionizing photons 

The equivalent width in H/3 is a quantity that monotonically 
decreases with time within a given HII region, or even within 
a given star formation episode of an HII galaxy. The number 
of ionizing photons decreases with time, decreasing the H/3 
intensity, and the continuum at the same wavelength rises. The 
equivalent width in H/3 has been u sed ex tensively as a nebular 
age in dicator by Stasin ska et al. Stasinska & Izotov 

i2003h . Terlevich et al. d2004 and D ottori & Bica Jl98lh . as 
was first suggested by Dottori ( 1981). 

The underlying non-ionizing population, however, blurs 
this picture by raising the continuum; the H/3 equivalent width 
anti-correlates with oxygen abundance (though with a huge 
scatter) (Terlevich et al. (2004)), definitely pointing towards 
a contamination of the continuum by previous populations. 



Furthermore, stellar popul ation studies (e.g., Wes tera e t al. 
d2004h : Kong et al. 120031) : Cid-Fernandes et al. (2003)) in- 
dicate that stellar formation events last for times of the order 
of at least 100 Myr, not 10 Myr (which is the photoionization 
timescale). A stellar formation event in an HII galaxy naturally 
accumulates an intermediate age population (between 10 and a 
few 100 Myr) that participates in the continuum budget. 

The two effects (i.e. (i) the accumulation of old popula- 
tions from previous bursts, and (ii) the accumulation of an in- 
termediate age population due to the present burst) have op- 
posite effects in terms of interpreting H/3 equivalent widths as 
star-forming event age indicators. The first one only introduces 
scatter in the relation, while the second one helps to decrease 
H/3 equivalent widths with increasing burst ages. This latter ef- 
fect happens, however, in a manner that may not be as simple 
to evaluate as can be naively done by simply comparing with 
photoionization by single aged stellar populations. 

Figure |6] shows a plot of AlogQ(H°) versus the equivalent 
width in H/3, as calculated in the previous subsection. As in 
Fig. |4] we can see a clustering of most data in AlogQ(H°) 
around a positive value of 0.8. A big scatter exists as well, 
and as we will explain in Sect. 3.3., we believe that we 
can explain it just from the uncertainties in the population 
synthesis process. However, in contrast to Fig. 0] we see a 
strong correlation of the clustered data with the equivalent 
width in H/3. A log Q(H°), which can be considered an indica- 
tor of the number of ionizing photons escaping the nebulae, 
increases with the decreasing equivalent width in H/3, i.e., time. 

This result may seem counter-intuitive, as the quantity 
of matter absorbing the ionizing photons remains the same 
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Fig. 6. Logarithmic ratio of the predicted number of ionizing 
photons to the number of recombinations as a function of the 
equivalent width of H/3. 

and the number of ionizing photons decreases with time: the 
optical thickness should increase with time. However, if we 
think of HII regions not as homogeneous media, but consider 
them in their complexity, with all their clumps, filaments, 
and substructures, we can understand the phenomenon, if we 
assume that the nebulae are optically thin in some directions, 
and thick in others. This can be described by a covering factor 
f c , characterizing the fraction of solid angle within which 
the nebulae are optically thick. The increase of the number 
of ionizing photons escaping the nebulae with time can thus 
be understood if the covering factor f c decreases with time . 
This has already been pointed out by Stasinska et al. (2001) 
and Stasinska & Izotov (2003), based on other arguments (the 
time sequence of the intensities of collisional lines in HII 
galaxies). The decrease of the covering factor f c with time can 
be understood if the structure of HII regions always becomes 
more complex with time, continuing to subfragment and 
thus to form more and more "holes" in their structure while 
they expand, always letting more ionizing photons escape. 
Tenorio-Tagle et al. J200 6) reach similar conclusions from the 
theoretical modelling of the hydrodynamical evolution of giant 
HII regions under the effect of both photoionization and stellar 
winds in a clumpy interstellar medium. 

Of course, the absorption of ionizing photons by dust would 
provide an alternative explanation; it would, however, require 
a mechanism that increases the absorption of ionizing photons 
by dust with time. As such a mechanism seems rather artificial 
to us, we prefer to stick to the leaking photons explanation. 

The biggest problem with leaking photons is that they are 
not confirmed by the direct observation of the Lyman contin- 



uum in local starburs t galax ies, at least in the quantities we pre- 
dict (Leitherer et al. (1995); Heckman et al. J2001I) : Deharveng 
et al. d200ll) : Bergvall et al. d2006l) ^. although is not the case at 
higher redshifts (Steidel (2001)). The total number of galaxies 
that have been observed in their Lyman continuum is, however, 
quite small, just a little over 10. These galaxies have been se- 
lected on the basis of the intensity of their starbursts. In the 
explanation we propose, the youngest starbursts, i.e., ones with 
the highest H/3 equivalent widths, which would be interpreted 
as the ones having the most intense starbursts, have the lowest 
fraction of escaping ionizing photons. We evaluate this fraction 
in these galaxies to be ^70%, i.e., AQ(H°) 3 on a linear scale; 
one has to take into account that our evaluations of AQ(H°) are 
certainly systematically over-evaluated by aperture effects by a 
factor that can be as high as 3 (see Sect. 2). Adding some 
systematic absorption by dust (which we, however, do not be- 
lieve can be much higher than 20%), we can easily lower the 
actual escaping fraction of ionizing photons in those galaxies 
to si 10%, which is what is believed to be the upper limit for 
escaping photons from Lyman continuum observations. This 
is, however, certainly not the case for the most evolved, low- 
est H/3 equivalent width galaxies in our sample, where AQ(H°) 
reaches values higher than 10. 

3.3. Accuracy of the population synthesis modelling 

The interpretation of the data presented above highly relies on 
our ability to evaluate the spectrum in the EUV from a popula- 
tion synthesis over the observed range between 4000-7000 A. 
Population synthesis is an ill-conditioned problem, subject to 
all kinds of degeneracies, and it is important to test its ability 
to recover what we want. 



3.3.1 . Comparison of the different libraries 

We first compared the results given by three different stel- 
lar population codes (or combinations), "BC03", "BC99", and 
"Starburst". The results commented on in Sect. I3.2I are overall 
consistent with the three different codes, but some differences 
do exist. It is noteworthy that, as far as the correlations with 
the H/3 equivalent widths are concerned, the best results are ob- 
tained with "Starburst" (Fig.|6}. "BC99" gives nearly as good 
correlations, and with fewer outliers. Finally, results obtained 
with "BC03" do not give any improvement; on the contrary, 
the correlation presented, though definitely existing, is not as 
strong as that with the two other libraries, and the number of 
outliers is higher. 

Though this cannot be considered definite evidence, since 
effects on the number of ionizing photons is what we want to 
assess, this could be indicative of deficiencies of the BC03 
stellar spectra library in the EUV. The BC03 library is made 
up of the extension of th e emp irical STELIB library in the 
visible (Le Borgne et al. (2003)) to other wavelengths, with 
empirical spectra when available, and using theoretical spectra 
when empirical spectra were not available. The connection 
of the different spectra of different origins has been based on 
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evaluated physical parameters, or rather, parameter classes, 
introducing some inconsistency in the process, which is not 
present in entirely theoretical libraries; one must remember 
that the EUV part of the spectra is theoretical in all libraries. 



3.3.2. Monte-Carlo simulations 

As a matter of fact, in our present problem, we are not so much 
interested in the underlying population parameters, as in our 
ability to recover a given spectrum from observations over a 
limited wavelength range. 

To do this, we took the synthesised spectra over the whole 
91 A to 160 fim, corresponding to the population mixtures de- 
rived from the population synthesis over the observed range 
(4000-7000 A). We added Gaussian noise, such that the S/N 
ratio varied from 2 to 40, corresponding to the observed S/N 
ratio distribution (see Kehrig et al. (|2004)), to the fitted syn- 
thetic spectrum (assumed to be our true prior), and recovered 
new synthesised spectra through the same population synthe- 
sis process over the 4000-7000 A range. We assumed the first 
spectra to give the true number of ionizing photons, giving the 
number of recombinations counted in H/3 under the hypothe- 
sis of the nebulae's optical thickness. From the latter spectra, 
we compute the number of ionizing photons by the population 
synthesis process, and we obtain A log Q(H°) from the logarith- 
mic difference of both numbers. To obtain reliable statistics, we 
performed 1620 such Monte-Carlo simulations for each of the 
three SSP libraries. 

Figures □andODrepresent A log Q(H°) as a function of oxy- 
gen abundance and Hfi equivalent width, respectively. If the 
population synthesis process were exact, the recovered num- 
ber of ionizing photons would be equal to the true number of 
ionizing photons, and in both cases A log Q(H°) would be con- 
stant, equal to 0. The only thing we measure here is the ability 
of the population synthesis process to recover the number of 
ionizing photons. As can be seen in Figs. and [8] the disper- 
sion obtained is equivalent to the dispersion observed in Figs.0] 
and|6j all the scatter in AlogQ(H°) can simply be explained 
by the uncertainty introduced by the population synthesis pro- 
cess. However, as can also be seen in Figs. [71 and [S] the popula- 
tion synthesis process does not introduce any systematical bias: 
most values of A log Q(H°) lie around 0, and the distribution is 
fairly symmetric. The systematical clustering of the values of 
A log Q(H°) around 0.8 cannot be explained simply by an effect 
of the population synthesis process, and must be real. 

4. Summary 

For a sample of HII galaxies, we compared the number of pre- 
dicted ionizing photons from the population synthesis of their 
spectra in the visible to the number of recombinations counted 
in the observed Hfi. 

Differences between these two numbers can arise (i) be- 
cause the evaluation of the number of ionizing photons is erro- 
neous, either because the theoretical ionizing spectra are wrong 
or because the population synthesis algorithm is not able to 
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Fig. 7. Monte Carlo simulation of the logarithmic ratio of the 
predicted number of ionizing photons to the number of recom- 
binations as a function of the oxygen abundance from the gas 
for optically thick nebulae. 



evaluate the ionizing spectrum correctly, or (ii) because a frac- 
tion of the ionizing photons escape since the nebulae is opti- 
cally thin. 

We checked for systematic tendencies depending on metal- 
licity (as measured by the oxygen abundance from the gas) that 
would be indicative of systematic deviations of the theoretical 
spectra in the EUV, as far as the number of hydrogen ionizing 
photons is concerned. We checked this with three different pop- 
ulation synthesis libraries, and we did not find any systematic 
effects for our sample. 

Using Monte-Carlo simulations, we checked for the abil- 
ity of our population synthesis method to recover the actual 
number of ionizing photons. Though most cases provided great 
scatter, we did not detect any systematic trend. 

Finally, we examined the behaviour of the ratio of the num- 
ber of predicted ionizing photons to the number of recombi- 
nations as a function of Hfi equivalent width. Though affected 
by the presence of older populations, the equivalent width in 
H/3 can be considered a gross indicator of the nebular age, de- 
creasing as the HII nebula becomes older. We found a strong 
tendency for an increasing number of ionizing photons escape 
the nebulae as the nebulae get older. 

We interpret this as a consequence of the geometrical evo- 
lution of the HII regions, which should be optically thick to the 
ionizing photons in some directions and optically thin in oth- 
ers. We can understand it in terms of a covering factor f c of the 
ionized gas that should decrease with time. In other words, HII 
nebulae subfragment more and more as they expand, always 
letting more ionizing photons escape. 




10 



F. Cuisinier et al.: On the geometrical evolution of the ionized gas in HII galaxies 



1 1 1 1 1 1 1 1 1 1 1 

"BC99"- 

_ " * ' X " 

X K s , * 

|V *, « , ; 


i i i i i i i i i i i 

"Starburst"- 

>f "3?" It- * * w * ^ 


1 1 1 1 1 1 1 1 1 1 1 
"BC03"- 


- * " - 


■ * x< ' xx 

K x 4 ' - 

" 


_Jx fc'Sx* '^x^ , 
X* x * 

_ x - if ** x »* 



1.5 2 1.5 2 1.5 2 



log(EW(H,)) log(EW(H„)) log(EW(H,)) 
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binations as a function of the equivalent width in H/? for opti- 
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The absorption of ionizing photons by dust would provide 
an alternative explanation, but we do not favour it, as it would 
require dust to absorb an always increasing fraction of ionizing 
photons during the evolution of HII regions. 

Acknowledgements. We would like to thank FAPERJ, CNPq, and the 
Swiss National Foundation for financial support. F. C. wishes to thank 
the Astronomisches Institut der Universitat Basel for hospitality dur- 
ing a stay during which most of this paper was written. We wish to 
thank the anonymous referee for very helpful comments that helped 
to clarify and expand the discussion. 

References 

Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 
Bergvall, N., Zackrisson E., Andersson, B.-G., et al., 2006 
A&A,448,513 

Bland-Hawthorn, J., & Maloney, P. R. 1999, ApJ, 510, L33 
Bressan, A., Fagotto, E, Bertelli, G, & Chiosi, C. 1993, 
A&AS, 100, 647 

Bruzual A., G, & Chariot, S. 1993, ApJ, 405, 538 
Bruzual A., G, & Chariot, S. 2000, Galaxy isochrone spectral 
synthesis evolution library (private communication) 
Bruzual A., G., & Chariot, S. 2003, MNRAS, 344, 1000 
Calzetti D., Kinney A. L., & Storchi-Bergmann T. 1994, ApJ, 
429, 582 

Castellanos, M., Dias, A. I., & Tenorio-Tagle, G. 2002, ApJ, 
565, 79 

Cervino, M, & Valls-Gabaud, D. 2003, MNRAS, 338, 481 
Charbonnel, C, Maeder, A., Schaller, G, et al. 1993, A&AS, 
101,415 



Chariot, S., & Bruzual A., G. 1991, ApJ, 367, 126 

Cid Fernandes, R., Leao, J. R. S., & Lacerda, R. R. 2003, 

MNRAS, 340, 29 

Cuisinier, F, Telles, E., & Kehrig, C. 2006, in preparation 
Deharveng, J.-M., Buat, V., Lebrun, V., et al. 2001 , A&A, 375, 
805 

Dottori, H. 1981, Ap&SS, 80, 267 

Dottori, H., & Bica, E. 1981, A&A, 102, 245 

Fagotto, F, Bressan, A., Bertelli, G., & Chiosi, C. 1994a, 

A&AS, 104, 365 

Fagotto, F, Bressan, A., Bertelli, G., & Chiosi, C. 1994b, 
A&AS, 105, 39 

Ferland, G. 1986, PASP, 98, 549 

Fluks, M. A., Plez, B., The, P. S., de Winter, D., Westerlund, 
B. E., & Steenman, H. C. 1994, A&AS, 105, 311 
Girardi, L., Bressan, A., Chiosi, C, Bertelli, G, & Nasi, E. 
1996, A&AS, 117, 113 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, 

A&AS, 141, 371 ("Padova2000" isochrones) 

Heckman, T, Sembach, K., Meurer, G., et al. 2001, ApJ, 558, 

56 

Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407 
Hirashita, H, Buat, V., & Inoue, A. K. 2003, A&A, 410, 403 
Inoue, A.K., 2001, AJ, 122, 1788 

Inoue, A.K., Hirashita, H., & Kamaya, H. 2001, AJ, 555, 613 
Kehrig, C, Telles, E., & Cuisinier, F. 2004, AJ, 128, 1141 
Kong, X., Chariot, S., Weiss, A., & Cheng, F. 2003, A&A, 
403, 877 

Le Borgne, J.-F, Bruzual A., G, Pello, R., et al. 2003, A&A, 
402, 433 

Leitherer, C, Ferguson, H. C, Heckman, T, & Lowenthal, J. 
D. 1995, ApJ, 454, L19 

Leitherer, C, Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 
123,3 

Lejeune, T, Cuisinier, F, & Buser, R. 1997, A&AS, 125, 229 
Lejeune, T, Cuisinier, F, & Buser, R. 1998, A&AS, 130, 65 
Meynet, G, Maeder, A., Schaller, G., Schaerer, D., & 
Charbonnel, C. 1994, A&AS, 103, 97 
Morisset, C. 2004, ApJ, 601, 558 

Mowlavi, N, Meynet, G, Maeder, A., et al. 1998, A&AS, 
128,471 

Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 2001, 
A&A, 375, 161 

Pequignot D. 1980, A&A, 83, 52 
Pickles, A. J. 1998, PASP, 1 10 863 

Sampson, L. 2003, Historia de Formacao Estelar em Galaxias 
HII, MSc. dissertation, Observatorio Nacional, Brazil 
Schaerer, D., Charbonnel, C, Meynet, G, Maeder, A., & 
Schaller, G. 1993a, A&AS, 102, 339 

Schaerer, D., Meynet, G., Maeder, A., & Schaller, G. 1993b, 
A&AS, 98, 523 

Schaerer, D., & de Koter, A. 1997, A&A, 322, 598 
Schaerer, D., & Vacca, W. D. 1998, ApJ, 497, 618 
Schaller, G., Schaerer, D., Meynet, G, & Maeder, A. 1992, 
A&AS, 96, 269 

Schaller G., Schaerer, D., Meynet, G, & Maeder, A. 1993, 
A&AS, 98, 523 



F. Cuisinier et al.: On the geometrical evolution of the ionized gas in HII galaxies 



11 



Schmutz, W., Leitherer, C, & Gruenwald, R. 1992, PASP, 
104, 1164 

Smith, L. J., Norris, R. P. R, & Crowther, P. A. 2002, MNRAS, 
337, 1039 

Stasinska, G. 1978, A&AS, 32, 429 

Stasinska, G., & Schaerer, D. 2001, A&A, 370, 1 

Stasinska, G., & Izotov, Y. 2003, A&A, 397, 51 

Steidel, C. C, Pettini, M., & Adelberger, K.G. 2001, ApJ, 546, 

665 

Tenorio-Tagle, G, Munoz-Tunon, C, Silich, S., Perez, E., & 
Telles, E., 2006, ApJ, in press 

Terlevich, R., Melnick, J., Masegosa, J., Moles, M., Copetti, 
M. V. E, 1991, A&AS, 91, 285 (T91) 

Terlevich, R., Silich, S., Rosa-Gonzales, D., & Terlevich, E. 
2004, MNRAS, 348, 1191 

Tylenda, R., & Stasinska, G. 1989, A&A, 217, 209 
Westera, P. 2001, The BaSeL 3.1 models: Metallicity cali- 
bration of a theoretical stellar spectral library and its appli- 
cation to chemo-dynamical galaxy models, PhD thesis, Univ. 
of Basel, 378 pp. 

Westera, P., Lejeune, T, Buser, R., Cuisinier, E, & Bruzual 
A., G. 2002, A&A, 381,524 

Westera, P., Cuisinier, E, Telles, E., & Kehrig, C. 2004, A&A, 
423, 133 

Worthey, G., Faber, S. M., Jesus Gonzalez, J., & Burstein, D. 

1994, ApJS, 94, 687 

Worthey, G. 1994, ApJS, 95, 107 

Zanstra, H. 1931, Zeitung fur Astrophysik, 2, 1 

Zurita, A., Rozas, M., & Beckman, J. E. 2000, A&A, 363, 9 



Appendix A: Hp equivalent widths and population 
parameters of individual spectra 



12 



F. Cuisinier et al.: On the geometrical evolution of the ionized gas in HII galaxies 




2 3 4 5 

log (A) [A] 



Fig. 2. Synthesised spectra over the whole wavelength range. The thin (red) lines represent the empirical spectrum (as in Fig.Q, 
the solid lines show the best-fitting spectra, and the dashed lines show the decompositions of the best-fitting spectra into young, 
intermediate, and old populations, respectively. The shaded regions represent the radiation shortward of the Lyman limit (912 
A), which is responsible for the ionization of the surrounding gas. Note that both axes are on a logarithmical scale, in contrast to 
Fig.H] 
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Table A.l. equivalent widths, A log Q(H°) parameters, mass ratios (relative to the total stellar mass), and ages of the individual 
sub-populations for each spectrum. 
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- 


0.9901 


CTS1013 


129.8 


0.679 


0.0029 


le+06 


0.0293 


5.088e+08 


0.9677 


CTS1016 


126.9 


0.952 


0.0083 


le+06 


0.0826 


2.026e+08 


0.9091 


CTS1018 


309.6 


0.859 


0.0161 


le+06 


0.0161 


5.088e+08 


0.9677 


CTS1019 


423 


-0.221 


0.0025 


3.981e+06 


0.0074 


5e+07 


0.9901 


CTS1028 


344.4 


0.588 


0.0227 


le+06 


0.0682 


5e+07 


0.9091 


CTS1033 


259.4 


0.409 


0.0029 


le+06 


0.0293 


2.026e+08 


0.9677 


CTS1034 


94.08 


0.859 


0.0029 


1.995e+06 


0.0293 


1.015e+08 


0.9677 


DDO060(NW) 


30.67 


-1.637 





- 


0.0909 


2.026e+08 


0.9091 


DDO060(SE) 


38.55 


-1.751 





- 


0.0909 


2.026e+08 


0.9091 


DDO070(E) 


225.3 


0.796 


0.0081 


le+06 


0.0242 


5.088e+08 


0.9677 


DDO070(W) 


210.9 


0.741 


0.0081 


le+06 


0.0242 


2.026e+08 


0.9677 


DDO075(CentSW) 


81.24 


0.831 


0.0161 


5.012e+06 


0.0161 


1.995e+07 


0.9677 


DD0155(NE) 


75.51 


-0.506 


0.0029 


6.918e+06 


0.088 


1.015e+08 


0.9091 


DD0155(SW) 


102.8 


1.037 


0.0293 


le+06 


0.8798 


5.088e+08 


0.0909 


ESO289IG037 


54.44 


-0.319 


0.0029 


7.943e+06 


0.0293 


1.015e+08 


0.9677 


ESO533G014 


32.28 


-1.648 





- 


0.0909 


2.026e+08 


0.9091 


FairalBO 


383.5 


0.556 


0.009 


le+06 


0.0009 


5.088e+08 


0.9901 


Haro24 


74.1 


-2.069 





- 


0.25 


2.026e+08 


0.75 


IC5154(N) 


40.17 


-0.162 


0.0083 


9.12e+06 


0.0826 


2.026e+08 


0.9091 


IC5154(S) 


53.58 


-0.495 


0.001 


6.918e+06 


0.0312 


1.015e+08 


0.9677 


MarseilleOl 


144.7 


-0.541 


0.0081 


9.12e+06 


0.0242 


5e+07 


0.9677 


Marseille68 


114.3 


0.893 


0.0025 


1.995e+06 


0.0074 


2.026e+08 


0.9901 


MBG0241 1-1457 


13.98 


-1.261 





- 


0.0909 


5.088e+08 


0.9091 


MBG21567-1645 


17.91 


-1.560 





- 


0.0323 


5e+07 


0.9677 


MBG22012-1550(E) 


62.18 


-0.409 


0.0081 


9.12e+06 


0.2419 


2.026e+08 


0.75 


Mrk36 


321.3 


0.727 


0.0227 


le+06 


0.0682 


5e+07 


0.9091 


Mrk710 


170.9 


0.937 


0.0081 


le+06 


0.0242 


5e+07 


0.9677 


Mrk710(CentNE) 


143.6 


-2.294 





- 


0.0909 


1.015e+08 


0.9091 


Mrk710(NE) 


185.1 


0.857 


0.0083 


le+06 


0.0826 


2.026e+08 


0.9091 


Mrk711 


135.9 


-0.065 


0.0003 


le+06 


0.0096 


1.015e+08 


0.9901 


Mrkl318 


323.9 


0.253 


0.0029 


1.995e+06 


0.0293 


5e+07 


0.9677 


NGC7323 


20.61 


0.366 





- 


0.0323 


le+07 


0.9677 


NGC7323(E) 


7.079 


0.337 


0.0003 


5.012e+06 


0.0096 


le+07 


0.9901 


NGC7323(W) 


13.28 


0.356 


0.0003 


5.012e+06 


0.0096 


le+07 


0.9901 


Tol0104-388(NW) 


267.6 


0.930 


0.0161 


1.995e+06 


0.0161 


5.088e+08 


0.9677 


Tol0104-388(SE) 


14.71 


0.218 







1 


5.088e+08 





Tol0117-414EW 


66.1 


-0.402 


0.0081 


9.12e+06 


0.2419 


5.088e+08 


0.75 


Tol0117-414NS(CentN) 


42.38 


-0.072 


0.0029 


6.918e+06 


0.088 


2.026e+08 


0.9091 


Tol0117-414NS(CentS) 


77.4 


0.907 


0.0029 


1.995e+06 


0.0293 


1.015e+08 


0.9677 


Tol0117-414NS(N) 


158.1 


-0.308 


0.0161 


7.943e+06 


0.0161 


1.995e+07 


0.9677 


Tol0226-390 


439 


0.324 


0.0081 


le+06 


0.0242 


1.995e+07 


0.9677 


Tol0306-405 


228.3 


0.749 


0.0161 


le+06 


0.0161 


5.088e+08 


0.9677 


Tol0341-407(E) 


167.3 


0.774 


0.0081 


le+06 


0.0242 


5e+07 


0.9677 


Tol0341-407(W) 


217.2 


-1.759 







0.0323 


1.995e+07 


0.9677 


Tol0440-381 


143.8 


1.194 


0.0323 


1.995e+06 







0.9677 


Tol0528-383(E) 


102.9 


-0.438 


0.0081 


7.943e+06 


0.2419 


5.088e+08 


0.75 


Tol0528-383(W) 


160.4 


-0.679 


0.0293 


9.12e+06 


0.8798 


5.088e+08 


0.0909 
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Table A.l. equivalent widths, A log Q(H°) parameters, mass ratios (relative to the total stellar mass), and ages of the individual 
sub-populations for each spectrum (continued). 



Spectra 


EW(H /J ) 


A log Q(H U ) 


My/M to , 


Age y 


Mi/M, B , 






M„/M tol 


Tol0610-387 


29.41 


-2.003 





- 


0.25 


2 


,026e+08 


0.75 


Tol0645-376 


122 


0.696 


0.0029 


1.995e+06 


0.0293 




5e+07 


0.9677 


Tol0957-278(NW) 


267 


0.662 


0.0081 


1.995e+06 


0.0242 




5e+07 


0.9677 


Tol0957-278(SE) 


186.5 


1.016 


0.0161 


le+06 


0.0161 


2 


.026e+08 


0.9677 


Toll004-296(NW) 


330.9 


0.683 


0.0074 


le+06 


0.0025 


2 


.026e+08 


0.9901 


Toll004-296(SE) 


306.3 


0.669 


0.0081 


1.995e+06 


0.0242 


1 


.015e+08 


0.9677 


Toll025-285 


53.61 


0.979 


0.0029 


1.995e+06 


0.088 


2 


.026e+08 


0.9091 


Toll 147-283 


168.9 


-0.349 


0.0025 


6.026e+06 


0.0074 


1 


.995e+07 


0.9901 


Toll 345-420 


226.8 


0.460 


0.0029 


le+06 


0.0293 


2 


.026e+08 


0.9677 


Toll457-262E 


413.2 


0.637 


0.0099 


le+06 







- 


0.9901 


Toll457-262W(Cent) 


146.3 


0.645 


0.0029 


1.995e+06 


0.0293 


1 


.015e+08 


0.9677 


Toll457-262W(CentE) 


180.9 


0.887 


0.0083 


le+06 


0.0826 


5 


.088e+08 


0.9091 


Toll457-262W(CentW) 


236.7 


0.717 


0.0081 


le+06 


0.0242 


1 


.995e+07 


0.9677 


Toll457-262W(E) 


204.5 


0.716 


0.0081 


le+06 


0.2419 


5 


.088e+08 


0.75 


Toll457-262W(W) 


215.4 


0.993 


0.0161 


le+06 


0.0161 


5 


.088e+08 


0.9677 


Toll 924-4 16(E) 


541.3 


0.642 


0.0682 


1.995e+06 


0.0227 




le+07 


0.9091 


Toll 924-4 16(W) 


497.2 


0.606 


0.0161 


le+06 


0.0161 




5e+07 


0.9677 


Toll 937-423 


31.84 


-0.375 


0.0029 


9.12e+06 


0.088 


2 


.026e+08 


0.9091 


UM69(Cent) 


46.25 


0.186 


0.0227 


9.12e+06 


0.2273 


5 


.088e+08 


0.75 


UM69(W) 


91.47 


-1.153 





- 


0.5 


1 


,995e+07 


0.5 


UM137(CentE) 


18.56 


-1.559 





- 


0.5 


5 


.088e+08 


0.5 


UM137(CentW) 


24.19 


-1.793 





- 


0.5 


5 


.088e+08 


0.5 


UM137(W) 


19.76 


1.309 


0.0029 


le+06 


0.088 


5 


.088e+08 


0.9091 


UM160(Cent) 


137 


0.409 


0.0161 


3.981e+06 


0.4839 


2 


.026e+08 


0.5 


UM 160(E) 


242.4 


-0.053 


0.0081 


5.012e+06 


0.0242 




le+07 


0.9677 


UM160(W) 


150.6 


0.277 


0.0083 


5.012e+06 


0.0826 


5 


.088e+08 


0.9091 


UM191 


37.62 


0.056 


0.0083 


9.12e+06 


0.0826 


5 


.088e+08 


0.9091 


UM238 


135.8 


-0.407 


0.0025 


6.918e+06 


0.0074 




le+07 


0.9901 


UM307 


118.1 


0.011 


0.0081 


6.026e+06 


0.2419 


5 


.088e+08 


0.75 


UM323 


113.6 


0.659 


1 


5.012e+06 







- 





UM395 


35.77 


0.154 


0.0081 


7.943e+06 


0.2419 


5 


.088e+08 


0.75 


UM396 


501.3 


0.484 


0.0161 


le+06 


0.0161 


1 


.015e+08 


0.9677 


UM408 


156.1 


0.942 


0.0081 


1.995e+06 


0.0242 


1 


.015e+08 


0.9677 


UM439(Cent) 


110.2 


0.979 


0.0083 


1.995e+06 


0.0826 


2 


.026e+08 


0.9091 


UM439(NW) 


223.3 


0.759 


0.0081 


1.995e+06 


0.0242 




5e+07 


0.9677 


UM448(extensionNE-SW) 


242.6 


0.400 


0.0029 


1.995e+06 


0.0293 




5e+07 


0.9677 


UM455(NW) 


190.2 


0.554 


0.0029 


1.995e+06 


0.0293 


2 


.026e+08 


0.9677 


UM456(Cent) 


239.9 


0.721 


0.0227 


le+06 


0.2273 


1 


.015e+08 


0.75 


UM456(NE) 


403.2 


0.720 


0.0682 


le+06 


0.6818 


5 


.088e+08 


0.25 


UM456(SW) 


215.5 


0.840 


0.0227 


le+06 


0.0682 


1 


,995e+07 


0.9091 


UM461(E) 


854.2 


0.448 


0.0323 


le+06 









0.9677 


UM461(W) 


638.4 


0.598 


0.0682 


le+06 


0.0227 


5 


.088e+08 


0.9091 


UM462(NE) 


301.5 


0.785 


0.0455 


le+06 


0.0455 


1 


.995e+07 


0.9091 


UM462(SW) 


382.4 


-0.158 


0.125 


5.012e+06 


0.375 




le+07 


0.5 


UM499(E) 


145.4 


-0.821 


0.0029 


9.12e+06 


0.0293 


2 


.026e+08 


0.9677 


UM499(W) 


58.06 


-0.529 


0.0029 


9.12e+06 


0.088 


2 


.026e+08 


0.9091 


UM533(E) 


138.4 


-0.807 


0.0161 


7.943e+06 


0.4839 


5 


.088e+08 


0.5 


UM559(Cent) 


33.93 


0.236 


0.0081 


7.943e+06 


0.2419 


2 


.026e+08 


0.75 


UM598W(CentNE) 


20.3 


-1.289 







0.0323 


2 


.026e+08 


0.9677 


UM598W(CentSW) 


46.34 


1.142 


0.0025 


1.995e+06 


0.0074 


2 


.026e+08 


0.9901 


UM598W(NE) 


58.18 


0.595 


0.001 


le+06 


0.0312 




5e+07 


0.9677 


UM598W(SW) 


112.8 


1.029 


0.0081 


le+06 


0.2419 


5 


.088e+08 


0.75 



